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(57) ABSTRACT 

Disclosed is a method for inter-antenna baseline determina- 
tion that uses an antenna configuration comprising a pair of 
relatively closely spaced antennas and other pairs of distant 
antennas. The closely spaced pair provides a short baseline 
having an integer ambiguity that may be searched exhaus- 
tively to identify the correct set of integers. This baseline is 
then used as a priori information to aid the determination of 
longer baselines that, once determined, may be used for 
accurate run time attitude determination. 

37 Claims, 13 Drawing Sheets 
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METHOD AND APPARATUS FOR RELIABLE 
INTER-ANTENNA BASELINE 
DETERMINATION 

FEDERAL GOVERNMENT LICENSE RIGHTS 5 

The U.S. Government has a paid-up license in this inven- 
tion and the right in limited circumstances to require the 
patent owner to license others on reasonable terms as 
provided for by the terms of NAS 1-20644 awarded by 
NASA. 10 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates generally to automated methods of 
attitude determination and particularly to methods employ- 15 
ing antenna arrays detecting radio waves. 

2. Background 

Pairs of antennas capable of sensing radio source carrier 
wave phases may be used to measure the attitude of an 
object. If the locations of the radio sources are known along 20 
with the number of cycles between antenna pairs, then 
detection of multiple radio sources reveals the attitude of an 
antenna array through ordinary geometry. An example of a 
device capable of producing radio source carrier waves is a 
satellite. 25 

Because the number and location of GPS satellites is 
known, they are well-adapted for use in making attitude 
determinations in the out-of-doors. However one is ordi- 
narily only able to sense the instantaneous phase angle of the 
radio source carrier wave. Thus for an antenna pair, with 30 
each antenna reporting a detected phase, there is an ambigu- 
ous number of cycles between the antennas. The ambiguity 
must be resolved in order for attitude to be determined from 
the radio sources. 

35 

The geometry of the antenna configuration limits the 
scope of the ambiguity in practice. However, the desire for 
long antenna baselines to achieve accurate attitude measure- 
ments results in the integer ambiguity associated with such 
baselines being computationally impractical to resolve 
through the evaluation of every candidate integer set. Unre- 
movable noise in carrier phase measurements further com- 
plicates integer ambiguity resolution. 

Establishing an acceptable integer ambiguity resolution 
system is vital to the usability of any GPS attitude 45 
instrument, however. Correct integer ambiguity resolution is 
most vital at the initial startup of the instrument, because 
before any attitude solutions can be calculated, the integers 
must be determined for at least a minimum set of satellites 
and antenna pairs. 5Q 

There are two ways in which integer ambiguity resolution 
can fail. The first is that no set of integers which meets the 
criteria for acceptance is found. This type of failure is not 
disastrous since it is possible to try again with another epoch 
of data. 55 

The second type of failure is wholly unacceptable, and 
that is when a set of integers is selected which meets the 
acceptance criteria, but which is in fact an incorrect set. Such 
false integer sets would lead to attitude determinations that 
would be in error, and if undetected, could result in the 60 
reporting of hazardously misleading information to the user. 

Thus it is critical that above all else, no false set of 
integers ever be allowed to be used in calculating GPS 
attitude that is subsequently passed on to the user of the 
instrument. 65 

In order to facilitate the discussion of integer search 
reliability, it is useful to define several measures of the 
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quality of the results obtained. The first measure is the pass 
rate, or the percentage of measurement sets for which an 
integer set is determined that passes all quality checks and 
is subsequently used to produce an attitude solution. The 
second measure is the true initialization rate, defined as the 
percentage of solutions that both pass and are correct. The 
final measure is the false initialization rate, defined as the 
percentage of passing solutions that are incorrect. 

It is relatively straightforward, as mentioned previously, 
to use a combination of quality checks to achieve a high pass 
rate and high true initialization rate. The very difficult matter 
is constructing a technique that preserves high pass rates and 
true initialization rates while suppressing false initializations 
entirely, or reducing the false initialization rate to zero. 

A related approach to integer ambiguity resolution can be 
presented as an example. There, for a particular antenna 
array consisting of three antennas and therefore three 
baselines, a double-difference integer ambiguity resolution 
scheme was employed that utilized appropriate limits on 
baseline magnitudes and inter-baseline angles. The candi- 
date integer sets were rank ordered according to the mag- 
nitude of the residual vector from the least-squares solution. 
Extensive efforts were made to select the best possible 
parametric values for the baseline magnitude and angle 
limits. As would be expected, the tighter the limits imposed, 
the lower the pass rates become. However, false initializa- 
tion rates, while reduced to less than 1%, were never entirely 
eliminated by this approach. 

The answer as to why such an approach cannot com- 
pletely eliminate false solutions is that in some particular 
situations of satellite - baseline geometry, so long as there is 
noise present in the carrier phase measurements, an incorrect 
set of integers will appear by all measures to represent a 
superior solution to the correct set of integers. This unfor- 
tunate circumstance is generally true whenever there is any 
integer ambiguity; a false initialization rate of zero cannot be 
guaranteed through the use of rank ordering and enforce- 
ment of solution quality criteria alone. 

There are at least two general classes of approaches to 
addressing this dilemma. 

The first approach involves the use of multiple epochs of 
data to determine the correct integer set. While it might be 
impossible to identify a false set of integers as false and 
select the true integer set given a single epoch of data, it will 
eventually become clear that the true set is superior to the 
false set as more epochs of data are gathered, particularly as 
the satellite geometries and/or baseline orientations change. 

A second method that can be used to aid the resolution of 
ambiguity is a priori knowledge of attitude. In fact, if the 
initial attitude is known precisely, then the integers can be 
directly calculated for each baseline and no integer ambi- 
guity exists. This method can be useful if there are auxiliary 
sensors that can be used to initialize the GPS attitude, or 
after reacquisition of GPS signals following an outage 
during which attitude has been updated using inertial sen- 
sors. However, since a priori knowledge of attitude will 
never be exact and since there will always be noise in carrier 
phase measurements, this approach by itself cannot be 
guaranteed to work flawlessly. 

Therefore a method and apparatus for reliable inter- 
antenna baseline determination is needed. 

SUMMARY OF THE INVENTION 

In light of this need, one aspect of the invention provides 
a method for determining inter-antenna baselines from radio 
wave phase measurements. This method identifies a first, 
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second and third antenna; determines a first baseline 
between the first and second antenna using what is termed a 
“backward” integer search procedure; and determines a 
second baseline between the second and third antennas using 
what is termed a “forward” integer search procedure respon- 5 
sive to the first baseline. 

Another aspect of the invention provides a method for 
selecting a baseline from a plurality of baseline candidates. 
This method ranks the baseline candidates and identifies a 
highest and next highest ranked baseline candidate. A dif- 10 
ference between the highest and next highest ranked candi- 
date is then computed. If the difference is large enough 
(indicating a uniquely best baseline candidate) the highest 
ranked is selected as baseline. 

Yet another aspect of the invention provides a method for 15 
determining an inter-antenna baseline from estimates of two 
or more orientation angles of the baseline. This method 
determines a limited geometric region from the estimates of 
the two or more orientation angles, identifies a plurality of 
baselines candidate vectors in the geometric region corre- 20 
sponding to integer sets that resolve an integer ambiguity 
associated with the baseline, and selects the baseline from 
among the vectors in the limited geometric region. 

Yet another aspect of the invention is an antenna configu- 
ration comprising a first, second, and third antenna. The first 25 
and second antenna are separated by a distance less than or 
equal to 1.75 wavelengths of a predetermined wave. The 
second and third antenna are separated by a distance greater 
than or equal to 3.0 wavelengths of the predetermined wave. 

BRIEF DESCRIPTION OF THE DRAWINGS 30 

FIG. 1 depicts a 'microprocessor based system’ in accor- 
dance with an illustrative embodiment; 

FIG. 2 depicts an 'antenna configuration’ in accordance 
with an illustrative embodiment; 35 

FIG. 3 A depicts an 'integer ambiguity associated with an 
inter-antenna baseline’ in accordance with an illustrative 
embodiment; 

FIG. 3B depicts an antenna-receiver unit in accordance 
with an illustrative embodiment; 40 

FIG. 4 depicts an 'inter- antenna baseline determination’ 
process in accordance with an illustrative embodiment; 

FIG. 5A depicts a 'first baseline determination’ process in 
accordance with an illustrative embodiment; 

FIG. 5B depicts a 'first baseline candidate integer set 45 
selection’ procedure in accordance with an illustrative 
embodiment; 

FIG. 5C depicts a 'first baseline correct candidate integer 
set selection’ procedure in accordance with an illustrative 
embodiment; 

FIG. 6A depicts a 'second baseline determination’ process 
in accordance with an illustrative embodiment; 

FIG. 6B depicts a 'second baseline candidate integer set 
selection’ procedure in accordance with an illustrative 55 
embodiment; 

FIG. 6C depicts a 'second baseline correct candidate 
integer set selection’ procedure in accordance with an illus- 
trative embodiment; 

FIG. 7A depicts a 'third baseline determination’ process in 60 
accordance with an illustrative embodiment; 

FIG. 7B depicts a 'third baseline candidate integer set 
selection’ procedure in accordance with an illustrative 
embodiment; 

FIG. 7C depicts a 'third baseline correct candidate integer 65 
set selection’ procedure in accordance with an illustrative 
embodiment; and 
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FIG. 8 depicts a plurality of baseline candidate vectors in 
a limited geometric region in accordance with an illustrative 
embodiment. 

DETAILED DESCRIPTION OF THE 
INVENTION 

Notations and Nomenclature 

b,y Baseline vector representing position of antenna i 
relative to antenna j 
c Speed of propagation of carrier wave 
DD iy (n-l)xl matrix of sensed carrier phase double 
differences for antenna pair ij 
LOS^ line-of-sight unit vector pointing toward radio 
source K 

LOS nx3 matrix of line-of-sight unit vectors toward n 
radio sources 

dLOS (n-l)x3 matrix of differences between radio source 
line-of-sight unit vectors 

N £ j k Integer for antenna pair ij and radio source K 
N iy nxl matrix of single -difference integers for antenna 
pair ij 

N iy (n-l)xl matrix of double -difference integers for 
antenna pair ij 

L Known distance between antennas i and j 
P Pseudoinverse matrix [LOS r LOS] _1 LOS r 
P^ Pseudoinverse matrix [dLOS r dLOS] _1 dLOS r 
SD,y nxl matrix of carrier phase single differences for 
antenna pair ij 

SD,y s nxl matrix of sensed carrier phase single differences 
for antenna pair ij 

At Ci y Actual time of measurement at receiver j minus that 
at receiver i (“receiver clock bias”). 

At it y Receiver-antenna propagation delay for receiver- 
antenna unit i minus that for receiver- antenna unit j 
(“receiver line bias”) 

L Carrier wavelength 

v/ Measurement noise for antenna-receiver pair ij and 
radio source K 

v-- nxl matrix of measurement noise for antenna-receiver 
pair ij 

v iJo (n-l)xl matrix of double-difference measurement 
noise 

Carrier phase measured at antenna i for radio source 
K 

Carrier phase measured at antenna j for radio source 
K 

(4 Sensed phase difference for antenna-receiver 
pair ij and radio source K 

(+.— +A* Measured phase difference for antenna- 
receiver pair ij and radio source K 
Rate of change of range from observer to radio source 
K 

Equations 


<t., K -^=(LOS^b v )/X 

(1) 

SD, y =[LOS]b,/X 

(2) 

b,j=X[LOS _1 ]SD ( y 

(3) 

b, r XPSD, y 

(4) 

P=[LOS r LOS] -1 LOS r 

(5) 

e, r SD, r [LOS][P]SD, y 

(6) 
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(CPLD), a field programmable gate array (FPGA), or other 
device known in the computing arts to be able to perform the 
necessary processing. 

FIG. 2 illustrates an antenna configuration 200 in accor- 
dance with one embodiment of the invention. Depicted 
comprising the antenna configuration 200 are a first antenna 
Al, a second antenna A2, a third antenna A3, and a fourth 
antenna A4. Also shown are a first baseline b 32 , a second 
baseline b 43 , a third baseline b 21? and a set of remaining 
baselines {b 41? b 3:L ,b 42 }. In one embodiment the antennas are 
coplanar although one skilled in the art will quickly realize 
that this is not fundamental to the present invention which 
also encompasses non-coplanar antenna configurations 
without limitation. Table 1 lists the dimensions of the 
antenna configuration in one embodiment of the invention. 

TABLE 1 




Value 

Value 

(13) 

Dimension 

(cm) 

(GPS LI wavelengths) 

20 

1^321 

16.0 

0.84 

(14) 

1^421 

60.0 

3.15 

|b 2 i|= |b 41 | 

134.4 

7.06 


25 


30 


(15) 

(16) 

(17) 

(18) 

( 19 ) 

35 

( 20 ) 

(21) 40 

(22) 


45 


50 


Description of Figures 

FIG. 1 depicts a microprocessor based system 100 
capable of embodying aspects of the invention. Shown are 
a microprocessor 110, a storage 120 and an input/output 
system 130. The input/output system 130 is capable of 
receiving information from a set of antenna-receiver units 
140 such as the 'antenna-receiver 5 unit 330 discussed below 
in connection with FIG. 3B. Further illustrated is a media 
drive 170, such as a disk drive, CD-ROM drive, or the like. 
The media drive 170 may operate with a computer-usable 
storage medium 175 capable of embodying computer- 
readable code able to configure the microprocessor-based 55 
system to embody aspects of the invention. The input/output 
system 130 may also operate with a keyboard 150 or a 
display 160. The input/output system 130 may also com- 
municate through well-known data transmission ports such 
as a RS-232 connector 180 or a NMEA0183 connector 185. 

As one skilled in the art would recognize, there are 
numerous structures of programmed or programmable logic 
capable of being configured to embody aspects of the 
invention. These may include, without limitation, a general 
purpose computer having either one or plural processors, a 
special purpose computer, an application specific integrated 
circuit (ASIC), a complex programmable logic device 


60 


65 


In reference to FIG. 2 and Table 1, in one embodiment the 
dimension |b 32 | corresponds to the distance the second 
antenna A2 and the third antenna A3 are located apart Dl. 
Similarly, the dimension |b 42 | corresponds to the distance the 
third antenna A3 and the fourth antenna A4 are located apart 
D2. 

One skilled in the art will realize that the present invention 
is not limited to applications using the GPS LI carrier 
frequency. Rather, in other embodiments of the invention, it 
would be possible to use, for example, the LI (1575.42 
MHz) and L2 (1227.60 MHz) carrier frequencies either 
alone or in concert. One embodiment of the present inven- 
tion would use each of a plurality of GPS carrier 
frequencies, e.g. LI and L2. Such an embodiment would 
determine baseline candidates on the basis of the plurality of 
frequencies. Candidates that are favorably evaluated using 
various figures of merit, such as the magnitude of the 
residual of a least squares solution, for measurements taken 
at plural frequencies may be regarded as more reliable 
candidates for the actual baseline sought. Implementations 
of such a method would be able to provide greater robust- 
ness and accuracy in the determination of inter-antenna 
baselines and are within the scope of the present invention. 

An overview of an illustrative embodiment can be pro- 
vided by reference to FIG. 2. The illustrative embodiment 
uses what are termed herein a “backward 55 integer search 
procedure and a “forward 55 integer search procedure. The 
backward approach consists of a nested search through all 
the possible values of each integer, and results in each 
possible integer set being tested for its quality. It is a search 
carried out through the full range of integer space. By 
contrast, the forward approach consists of using a priori 
knowledge of the possible locations of the two antennas with 
respect to one another to calculate a subset of integer 
combinations for consideration. It is a search carried out 
through geometric space rather than integer space. 

The backward approach is appealing when there is no a 
priori knowledge of the relative locations of the antennas, or 
when the separation of the antennas is small enough to limit 
the number of possible integer combinations. The forward 
approach is generally superior when there is a good a priori 
estimate of the antennas 5 positions with respect to one 
another. Both methods have potential pitfalls that must be 
guarded against in practice. 
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A combination of these two approaches provides 
practical, robust integer ambiguity resolution in combination 
with particular antenna arrays. 

In the illustrative embodiment, the first baseline b 32 is 
determined using a backward integer search procedure. The 5 
fact that the first baseline b 32 is chosen to be relatively short 
makes this exhaustive search of the feasible candidate 
integer sets computationally practicable. If no candidate 
integer set results in a baseline passing quality tests, or if two 
comparable best candidate integer sets are found, another 10 
epoch of data is used until the first baseline b 32 is deter- 
mined. While the illustrative embodiment may, on occasion, 
use multiple epochs of data, one skilled in the art will 
recognize that it can successfully determine all needed 
baselines based on a single epoch of data. 15 

When the first baseline b 32 is determined estimates of the 
heading and roll of the antenna configuration are made based 
on the first baseline b 32 . Using this information in conjunc- 
tion with knowledge of the magnitude and relative orienta- 
tion of the second baseline b 43 , a forward integer search 20 
procedure is used determine the second baseline b 43 . A 
similar search is then conducted to determine the third 
baseline b 21 . After the first, second, and third baselines are 
determined, the integers for the set of remaining baselines 
{b 4 i,b 31 , b 42 } may be conveniently determined from simple 25 
mathematical sums, since the sum of integers around any 
geometrically closed circuit is zero. Once computed, least- 
squares solutions for the set of remaining baselines {b 41 , 
b 31 ,b 42 } are calculated. 

FIG. 3 A illustrates an 'integer ambiguity associated with 30 
an inter-antenna baseline 5 300. A radio wave source 305 is 
shown from which a radio source carrier wave 310 propa- 
gates down and to the right of the figure. Typically there is 
more than one radio wave source 305. In an illustrative 
embodiment of -the invention, the radio wave sources are 35 
Global Positioning System (GPS) satellites and the radio 
source carrier wave 310 is a GPS carrier wave. One skilled 
in the art would recognize that the designation radio wave 
source 305 and radio source carrier wave 310 do not limit 
the scope of the present invention to those portions of the 40 
electromagnetic spectrum associated with radio broadcasts 
and, rather, the present invention is applicable to other 
portions of the electromagnetic spectrum without limitation. 
One skilled in the art would further recognize that the 
present invention is not limited to use with orbiting radio 45 
sources. Rather the present invention may also be used with 
so called "pseudolites, 55 for example radio sources installed 
indoors, fixed outdoor sources, or mobile outdoor sources 
that are not in orbit. 

An 'antenna i 5 320 and an 'antenna j 5 340 are generically 50 
depicted as is a Tine of site 5 unit vector 315 ("LOS 55 ) 
pointing toward the radio wave source 305. 'Antenna i 5 320 
and 'antenna j 5 340 are arbitrarily chosen but may be, for 
example, antennas depicted in the antenna configuration 200 
of FIG. 2. They form an antenna pair ij. Depicted between 55 
'antenna i 5 320 and 'antenna j 5 340 is a baseline 325. The 
baseline 325 may be, for example, a baseline depicted in the 
antenna configuration 200 of FIG. 2. The baseline 325 
represents a vector parallel to, and of equal magnitude, as a 
line segment connecting ' antenna i 5 320 and ' antenna j 5 340. 60 
At the moment in time depicted in FIG. 3 A, 'antenna i 5 320 
is on a locus of phase 0.0 wavelength 350 and 'antenna j 5 
340 is on an a locus of phase 0.5 wavelength 360. 

Depicted in further detail in FIG. 3B is an 'antenna- 
receiver 5 unit 330 comprising an antenna 333, such as 65 
'antenna i 5 320 or 'antenna j 5 340. The antenna 333 is 
connected to a receiver 335 comprising a phase reporter 337, 
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such as a phase locked loop. The receiver 335 reports a 
phase 339. The 'antenna-receiver 5 unit 330 may contain a 
clock 338. The clock 338 may be common to multiple 
receiver-antenna units or multiple receiver-antenna units 
may each have their own non-common clock. While the 
present description is of a configuration in which for each 
antenna there is a receiver, one skilled in the art will 
recognize that others may be employed. For example, one 
embodiment would employ a receiver for antenna pairs, and 
this receiver would report phase differences between the 
antenna pair. Another embodiment would use a receiver to 
report the phase detected at multiple antennas, e.g. multi- 
plexing. When multiplexing would be used in this manner, 
a single receiver could be employed to report the phases that 
are present at each antenna in the antenna configuration. 

In the illustrative embodiment, the receivers attached to 
'antenna i 5 320 and 'antenna j 5 340 will report the fractional 
carrier phases of 0.0 and 0.5 wavelengths, respectively as the 
phases present for a particular time of measurement. 
However, the receivers can only detect instantaneous phase 
of the radio source carrier wave 310 at each antenna and 
cannot measure the whole number of cycles of the radio 
source carrier wave 310 that exist between 'antenna i 5 320 
and 'antenna j 5 340. Thus, there is an integer ambiguity 380 
associated with the baseline 325 comprising the ambiguous 
whole number of cycles of the radio source carrier wave 310 
which lie between 'antenna i 5 320 and 'antenna j 5 340. The 
integer ambiguity 380 thus creates an ambiguity in the exact 
number (whole plus fractional) of cycles of the radio source 
carrier wave 310 that lie between 'antenna i 5 320 and 
'antenna j 5 340. From the reported fractional phase differ- 
ence of 0.5 (0.5 minus 0.0) between 'antenna i 5 320 and 
'antenna j 5 340, the exact phase difference of the radio 
source carrier wave 310 could be 0.5+N where N is any 
integer value. A set of possible locations for antenna j 370 
corresponding to N=-l to 3 are depicted, wherein the actual 
location of 'antenna j 5 340 corresponds to N=2. 

As one skilled in the art would recognize, in the absence 
of measurement error or other bias, and assuming the 
distance from the radio wave source 305 is large in com- 
parison to the magnitude of the baseline 325, the exact 
number of cycles of the radio source carrier wave 310 may 
be used to reliably resolve the integer ambiguity 380 and 
thereafter determine the baseline 325. 

For example, assuming there are n radio wave sources 
such as radio wave source 305 sensed by 'antenna i 5 320 and 
'antenna j 5 340 without measurement error or other bias, the 
difference between the carrier phase measurements taken 
simultaneously at 'antenna i 5 320 and 'antenna j 5 340 is 
related to a simple dot product as follows: 

^-^(LOS^A (l) 

where 

(^/^Carrier phase measured at antenna i for radio 
source K (units of wavelengths) 

(|)^=Carrier phase measured at antenna j for radio 
source K 

b i;/ =Baseline vector representing position of antenna i 
relative to antenna j 

LOS^=Line-of-sight unit vector pointing toward radio 
source 

X= Carrier wavelength 

and K=l, . . . , n. 

Taking the measurements for n radio carrier wave sources, 
a set of n equations similar to Equation 1 results which may 
be rewritten as a matrix equation: 

SDy=[LOS]bjyA (2) 
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where 

SD iy =nxl matrix of carrier phase single differences for 
the antenna pair 

LOS=nx3 matrix of line-of -sight unit vectors toward 
the n radio sources. 

If n=3, this equation can be solved for the baseline 
vectorb iy : 

by=^'[LOS _1 ]SD I y (3) 

More generally, if n=^3, the least squares solution for b iy 
becomes 

b iy =XPSD iy (4) 

where the pseudoinverse matrix P is 

P=[LOS r LOS] -1 LOS r (5) 

and resulting from the least squares solution for the 
baseline vector b-j is a least squares residual: 

e r SD iy -[LOS][P]SD i7 (6) 

However, in practice measurement error and other bias 
typically are present. One type of bias is a line bias, At iiy . A 
clock bias At Ciy . can arise when the clock 338 is not common 
to each receiver. In one embodiment of the invention a 
non-common clock is used while another embodiment of the 
invention uses a common clock. When a common clock is 
used there is no clock bias or there is an identifiable constant 
or nearly constant clock bias. In order that one may under- 
stand and practice the invention using either a common or 
non-common clock, an embodiment with the clock bias will 
be discussed. One skilled in the art will recognize that 
assuming a constant clock bias is equivalent to illustrating 
an embodiment with a common clock for the purposes of 
describing, understanding and practicing the present inven- 
tion. 

A “corrected” or “sensed” phase difference is derived 
from the raw measured phase difference by accounting for a 
measurement error term that is related to clock differences 
but which cannot be eliminated by double differencing. 

If the measurements from a particular pair of antenna- 
receiver units are not taken at precisely the same time, an 
error source will appear in the measurement of phase single 
differences due to the fact that a constellation of radio 
sources (for example satellites) may be in motion. During 
the interval between measurements by two GPS receivers, 
for example, the GPS satellite constellation will have in fact 
moved enough to significantly affect the measured phase 
difference. 

Even between independent receivers, firmware features 
typically regulate the epoch to within 1 ms of the top of the 
GPS second, so clock bias between receivers is typically on 
the order of 1 ms or less. In 1 ms, a GPS satellite will move 
on the order of 3 m. While the resulting change in the line 
of sight vector to the satellite is of no practical consequence, 
the change in range to satellite has a nonnegligible effect on 
the phase difference measurement and therefore must be 
corrected for. 

The corrected or sensed phase difference is thus related to 
the raw measured phase difference by the expression 


* , K k r K (23) 

5 

where 

(4>^— 4> jr " fir )^= sense d phase difference for antenna- 
receiver pair ij and radio source K 
(<l>^-^)^ measured phase difference for antenna- 
10 receiver pair ij and radio source K 

r^=rate of change of range from observer to radio 
source K 

X= carrier wavelength 

At Cjy = ac tual time of measurement at receiver j minus 
15 that at receiver i (also known as “receiver clock 

bias”) 

While the satellite range rate r^ is typically much smaller 
than the speed of wave propagation c, it has the property of 
being satellite-specific. Therefore, double differencing will 
20 not eliminate its effect on the phase measurement; in practice 
the phase difference measurements for individual satellites 
must be corrected by this term before any further processing 
takes place. Correcting for this term is desirable even in the 
case of a common-clock architecture that results in an easily 
25 estimated but non-zero clock bias. 

The correction term 

T Alaj 

30 

can be adequately applied using even a coarse estimate of 
receiver clock bias. For example, for multiple independent 
GPS receivers, by taking the difference in GPS time solu- 
tions from independent receivers, one can arrive at a very 
coarse estimate of clock bias that is typically accurate to 
within about 300 ns. This uncertainty corresponds to about 
100 m of carrier wave propagation but a maximum of only 
about 1 mm of change in satellite range. Therefore, this 
4Q coarse clock bias estimate is adequate to make the needed 
correction to phase difference for satellite motion during the 
interval between receiver measurement epochs. 

There are numerous ways to estimate the satellite range 
rate for use in this correction term, including ephemeris 
45 calculations and using the carrier wave Doppler shift. In one 
embodiment, a difference in whole carrier phases between 
sequential measurements is used. 

Reconsidering Equation 1 with measurement error and 
bias included and separating terms to reflect the integer 
5Q ambiguity 380 the exact phase difference of the radio source 
carrier wave 310 between 'antenna T 320 and 'antenna j’ 340 
may be expressed as: 

- tf) = [<tf - $) s + {(At aj + AM] + Atf ■ - ^ (8) 

55 

where 

(t|) A- (|) ; A) 5 = sensed phase difference for antenna- 
receiver pair ij and radio source K 
60 C=speed of propagation of carrier wave 
X=carrier wavelength 

At Ciy =act ua l time of measurement at receiver j minus 
that at receiver i; the “clock bias” for the receiver ij 
pair 

65 At iiy =receiver-antenna propagation delay for receiver- 
antenna unit i minus that for receiver-antenna unit j; 
the “line bias” for the ij receiver- antenna pair 
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N^=integer for antenna pair ij and radio source K 
v iy A= measurement noise for antenna- receiver pair ij 
and radio source K 

Stacking all the measurements for n satellites yields 

5 

SDij = SD t j s + - (A t Cl j + A tuj) + Nij - v t j 


12 

where the fractional portion is calculated such that 
-0.5 < |/rac[(^f - <pj) s + £(A t C ij + Ar L y)]| < 0.5 

Then, 


where 10 

SD^=nxl matrix of sensed carrier phase single differ- 
ences for antenna pair ij 

N iy =nxl matrix of single -difference integers for 
antenna pair ij 

v iy =nxl matrix of measurement noise for antenna- 15 
receiver pair ij. 

One skilled in the art will recognize that while the 
foregoing has been illustrated in terms of single differences 
of sensed carrier phase, that double differences may be used 
without limitation. When using a common clock 20 
architecture, single differences are preferably used. 

One embodiment of the invention uses double differences. 
More specifically, subtracting the measurements corre- 
sponding to one radio wave source from other measurements ^ 
eliminates the clock bias and the line bias terms to yield: 


dd„.=dd„. +N„, 


VS ‘JD iJD 


( 10 ) 


where 30 

DDy s =(n-l)xl matrix of sensed carrier phase double 
differences for antenna pair 
NDy D =(n-l)xl matrix of double -difference integers for 
antenna pair ij 

v iy - D =(n-l)xl matrix of double-difference measurement 35 
noise. 

In either the case of single differences or double 
differences, a geometrically feasible set of integers may be 
established reflecting the feasible phase difference between 
an antenna pair separated by a known distance. 

In particular, if two antennas are separated from one 
another by a known distance, then for any particular carrier 
wave source emanating from any particular direction, the 
carrier phase difference for the pair of antennas is bounded 
by 45 


l^-^SL/X. (11) 

where 50 

L=Known distance between antennas i and j 
\= Carrier wavelength. 

Considering first the case of single differences, in the 
previous expression for phase difference, the term in brack- 
ets can be replaced with its fractional portion with no loss of 55 
generality, since the term N iy ^ can take on any integer value. 

(tf - *>j) = ), + C -(M aJ + bt w )] + N$ - (12) 

60 

If v,-, K is taken to be zero, bounds on N,-,^ can be 
established as follows: 

\4>f ~<Pj\ = |/rac(0? - <Pj) s + + touj) + A# | ^ L/X < 13 > g5 


|N/|SlA+0.5 (15) 

In the illustrative embodiment of the invention values of 
measurement noise v i} K are typically much less than a 
wavelength. In the worst practical case, the addition of noise 
can expand the range of possible integers to 

|N/|^lVX + 0.5 + |v/U flJC (16) 

Analogously for embodiments using double differences, 
the exact phase differences have representation: 

C^i^-^)- ( < l> i M - < l>y A/ ) =fra c[ (<bi K -tyj K )s- (tyi M -tyj M )s] + Nij K -Nij M -Vij K + 

v tj M (17) 


or 


DD/ M =fac[DD/ M 5 ]+N/ M -v/ M (18) 

where the double -difference quantities for antenna pair ij 
and radio source pair KM are defined as follows: 

dd 

M KM =fJ M 

L ^ij ~~ ij ij 

(19a-d) 

Now the integer search can be carried out by choosing 
candidate values of and proceeding in similar fashion 
to the method just outlined for the single-difference integer 
search. 

However, the range of integers that must be considered 
with a double-difference search technique is greater than that 
for the single-difference case. Analogous to single- 
difference integers, the following range of values for double- 
difference integers can be established: 


|N I y KA/ | 2L/A.+0.5 (20) 

In accordance with embodiments of the invention, a 
relatively short baseline greater than one-half of the wave- 
length 90 may be chosen, for instance L=0.8X. Assuming 
that |v iy ^|<0.2X, we have 

lN/1^0.8+0.5 (21) 

in the case of single differences or, in the case of double 
differences, 

N/ M |^l. 6+0.5 (22) 

Considering the implications of Equation 21, one skilled 
in the art will recognize that for this choice of L, there are 
three possible values for each integer, {-1,0,1}. For a single 
baseline then, the total number of possible integer combi- 
nations assuming there are n radio wave sources measure- 
ments is 3”, which for 8 radio wave sources is 6,561 
possibilities. In general, the number of possible integer sets 
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will be (r s ) n , where r s is the number of integers that fall 
within the range of possible single -difference integers. 

The resulting range for the double-difference integer 
search and the combinatorial consequences of a straightfor- 
ward search through this entire range can be summarized in 
an analogous way to those for the single -difference integer 
search. The number of possible double-difference integer 
combinations becomes (r^,)" -1 , where r D is the number of 
double-difference integers within the range. (The exponent 
is n-1 rather than n because there are one fewer double- 
difference integers than there are single-difference integers.) 
In accordance with an embodiment of the invention, if a 
particular satellite is chosen as the “master satellite”, say, 
satellite M, and the double differences are all derived by 
subtracting that particular measurement from the measure- 
ments corresponding to the other satellites, a convenient 
situation results. Because the double-difference integers are 
in fact each the difference of a particular satellite’s single- 
difference integer and the master satellite’s single -difference 
integer, all members of the set of double difference integers 
for a particular antenna pair must lie within a particular 
range of one another, that range being the range of possible 
master satellite single -difference integers. 

Using this fact, the integer search procedure can consist of 
r 5 subsearches corresponding to r 5 values of the master 
satellite’s single -difference integer. Each of these sub- 
searches comprises a search through (r nm)”' 1 single - 
difference integer combinations for the remaining n-1 sat- 
ellite measurements. The “non-master” integer range r^^ is 
the number of integers that must be included to ensure that 
all possible double -difference integer sets are examined, 
given the range of values of the single-difference master 
integer. The required value of r^^ is the same as r s for all 
antenna separations except those between x.25 and x.50 
wavelengths (x being any integral value), where a transition 
region exists. Using this approach, the number of double- 
difference integer combinations that must be searched 
changes from (r^,)” -1 to ^(r^^)” -1 , which is equal to (r 5 ) n 
for most antenna separations. 

For an embodiment:)f the invention sensing carrier 
sources from 8 satellites, the integer range and number of 
potentially possible candidate integer sets for single - 
difference and double-difference integers is summarized 
below through a range of inter-antenna separations L 


Number Number 

of Integer of Integer 
Set Set 

Candidates Candidates 
for for 

n = 8 n = 8 

Single- Double- (Single- (Double- 

differ- differ- Differ- Differ- 
ence ence ence ence 

Integer Integer Search) Search) 

L A Range Range r s r NM (r s ) n 

<0.25 0 Oil 1 1 

0.25- 0 -1, ... 1 1 3 1 2,187 

0.50 

0.50- -1, ... 1 -1, ... 1 3 3 6,561 6,561 

0.75 

0.75- -1, ... 1 -2, ... 2 3 3 6,561 6,561 

1.00 

1.00- -1, ... 1 -2, ... 2 3 3 6,561 6,561 

1.25 

1.25- -1, ... 1 -3, ... 3 3 5 6,561 234,375 

1.50 


-continued 








Number 

Number 

5 






of Integer 
Set 

of Integer 
Set 







Candidates 

Candidates 







for 

for 







n = 8 

n = 8 



Single- 

differ- 

Double- 

differ- 



(Single- 

Differ- 

(Double- 

Differ- 

10 


ence 

ence 



ence 

ence 


L/A 

Integer 

Range 

Integer 

Range 

r s 

r NM 

Search) 

( r s) n 

Search) 

r s( r NM) n 1 


1.50- 

-2, ... 2 

-3, ... 3 

5 

5 

390,625 

390,625 


1.75 
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1.75- 

-2, ... 2 

-4, ... 4 

5 

5 

390,625 

390,625 

2.00 








2.00- 

-2, ... 2 

-4, ... 4 

5 

5 

390,625 

390,625 


2.25 








2.25- 

-2, ... 2 

-5, ... 5 

5 

7 

390,625 

4,117,715 


2.50 







20 

2.50- 

2.75 

-3, ... 3 

-5, ... 5 

7 

7 

5,764,801 

5,764,801 


2.75- 

-3, ... 3 

-6, ... 6 

7 

7 

5,764,801 

5,764,801 


3.00 








3.00- 

-3, ... 3 

-6, ... 6 

7 

7 

5,764,801 

5,764,801 


3.25 







25 

3.25- 

3.50 

-3, ... 3 

-7, ... 7 

7 

9 

5,764,801 

33,480,783 


3.50- 

-4, ... 4 

-7, ... 7 

9 

9 

43,046,721 

43,046,721 


3.75 








3.75- 

-4, ... 4 

-8, ... 8 

9 

9 

43,046,721 

43,046,721 


4.00 








30 

FIG. 4 illustrates an ‘inter- antenna baseline determina- 
tion’ process 400. In accordance with an embodiment of the 
invention, the process 400 operates with the antenna con- 
figuration 200 and antenna-receiver units 330 previously 
35 discussed. The ‘inter- antenna baseline determination’ pro- 
cess 400 initiates at a ‘start’ terminal 405 and continues to 
receive an ‘epoch of data’ input block 410. Measured phases 
are reported from the antenna-receiver units, such as the 
‘antenna-receiver’ unit 330. Sensed phases are generated 
from the measured phases through Equation 23 and the 
40 sensed phases are reported in the ‘epoch of data’ input block 
410. Processing continues to a ‘backward search for first 
baseline’ procedure 420 that performs a backward search 
procedure to determine a first baseline as will be discussed 
in greater detail in connection with FIG. 5. If the ‘backward 
45 search for first baseline’ procedure 420 reports a failure to 
determine the first baseline, a ‘first baseline identified’ 
decision procedure 425 exits along its ‘no’ branch and 
processing returns to receive another ‘epoch of data’ input 
block 410. If the ‘backward search for first baseline’ proce- 
50 dure 420 successfully determines the first baseline, the ‘first 
baseline identified’ decision procedure 425 exits through its 
‘yes’ branch and processing continues to a ‘forward search 
for second baseline’ procedure 430. As will be discussed in 
greater detail in connection with FIG. 6, the ‘forward search 
55 for second baseline’ procedure 430 determines a second 
baseline or reports a failure condition. If the failure condi- 
tion is reported a ‘second baseline identified’ decision pro- 
cedure 435 exits through its ‘no’ branch and processing 
returns to await another ‘epoch of data’ input block 410. If 
60 the second baseline is determined, the ‘second baseline 
identified’ decision procedure 435 exits through its ‘yes’ 
branch and processing continues to a ‘forward search for 
third baseline’ procedure 440. As will be discussed in greater 
detail in connection with FIG. 7, the ‘forward search for 
65 third baseline’ procedure 440 determines a third baseline or 
reports a failure condition. If the failure condition is reported 
a ‘third baseline identified’ decision procedure 445 exits 
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through its ‘no’ branch and processing returns to await 
another 'epoch of data’ input block 410. If the third baseline 
is determined, the 'third baseline identified’ decision proce- 
dure 445 exits through its 'yes’ branch and processing 
continues to a 'remaining baseline determination’ procedure 5 
450 that determines the remaining baselines in the antenna 
configuration 200. The first, second, and third baselines 
having been determined, the integers that resolve their 
associated integer ambiguities are known and are used to 
compute the remaining baselines using ordinary geometry 10 
well known in the art and the least squares residual for the 
remaining baselines is calculated. 

A 'baseline configuration check’ procedure 460 next 
checks the integrity of the set of baselines. One embodiment 
of the invention checks whether the maximum magnitude of 15 
residual vector elements is less than 0.22 GPS LI wave- 
lengths and whether the mean value of residual vector 
elements for all baselines is less than 0.08 GPS LI wave- 
lengths. If these checks are not satisfied, the process 400 
begins on a new epoch of data 410. Otherwise, a 'configu- 20 
ration limits’ decision procedure 470 begins. Embodiments 
of the invention allow for limits to be placed on array 
attitude values. If this is desired, the 'configuration limits’ 
decision procedure 470 exits through its 'yes’ branch and a 
'level platform’ decision procedure 480 is entered. The 25 
'level platform’ decision procedure 480 tests whether the 
previously-determined baselines are consistent with a plat- 
form with which the antenna array is used being level within 
a predetermined range, for example +/-20 degrees. This test 
condition may be used as a final quality check for platforms 30 
which are nominally level, for example a ship. One skilled 
in the art will appreciate that a similar procedure could be 
used with a non-level fixed-attitude platform, for example a 
conventional landing gear aircraft. Conveniently, an incli- 
nometer could generate a nominal attitude for a non-level 35 
fixed attitude platform. 

If the previously-determined baselines are inconsistent 
with the platform being level within the predetermined 
range, the 'level platform’ decision procedure 480 exits 
through its 'no’ branch and a 'failure condition’ output block 40 
485 is reported for this epoch of data. The procedure then 
completes through an 'end’ terminal 490. 

One skilled in the art will recognize that a coarse estimate 
of two orientation angles is obtained through the determi- 
nation of the first baseline b 32 by the 'backward search for 45 
first baseline’ procedure 420. One skilled in the art will 
further recognize that a coarse estimate of two orientation 
angles may be otherwise obtained, for example from a solid 
state inclinometer and/or magnetometer, and the remaining 
functional elements depicted in FIG. 4 may be used with the 50 
coarse estimate so obtained. 

FIG. 5 A illustrates a 'first baseline determination’ process 
500 that is invoked by the 'backward search for first base- 
line’ procedure 420. Processing initiates at a 'start’ terminal 
510 and continues to a 'first baseline candidate integer set 55 
selection’ procedure 520 and a 'first baseline correct candi- 
date integer set identification’ procedure 530 that will be 
discussed in greater detail with respect to FIG. 5B and FIG. 

5C, respectively. Briefly here, all geometrically feasible 
candidate integers are selected and then searched for a 60 
uniquely best baseline candidate to identify as the integer set 
for the first baseline. The underlying rationale for this is that 
if there are two candidate baselines that are insignificantly 
different, both should be rejected as the correctness of 
neither can be assured. Processing continues to a 'successful 65 
identification’ decision procedure 585. If the 'first baseline 
correct candidate integer set identification’ procedure 530 
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has identified a uniquely best baseline candidate, the 'suc- 
cessful identification’ decision procedure 585 exits through 
its 'yes’ branch and a 'first baseline and candidate integer 
set’ output block 590 returns the identified first baseline and 
its associated candidate integer set. Otherwise, the 'success- 
ful identification’ decision procedure 585 exits through its 
'no’ branch and a 'failure condition’ output block 595 
reports the identification failure to the calling process 400 
that will cause the 'first baseline identified’ decision proce- 
dure 425 to exit along its 'no’ branch and receive the 'epoch 
of data’ input block 410. Processing completes through an 
'end’ terminal 599. 

FIG. 5B illustrates the 'first baseline candidate integer set 
selection’ procedure 520 invoked by the 'first baseline 
determination’ process 500. Processing initiates at a 'start’ 
terminal 521 and continues to a 'feasible integer sets deter- 
mination’ procedure 523 that determines candidate integer 
sets that are geometrically feasible given the known mag- 
nitude of the first baseline to be determined. Recall that the 
first baseline b 32 is determined using the backward search 
procedure. This means that a search though all of the 
possible values of each integer of geometrically feasible 
integer sets will be conducted. Each of the integer sets will 
be tested for the quality of the solution they would provide 
for the first baseline b 32 . Thus, all feasible integer sets need 
to be selected for further processing. An 'all feasible integer 
sets selection’ procedure 525 performs this function and 
processing completes through an 'end’ terminal 527. 

FIG. 5C depicts the 'first baseline correct candidate 
integer set selection’ procedure 530 invoked by the 'first 
baseline determination’ process 500. Processing initiates at 
a 'start’ terminal 531 and continues to a 'least squares 
computation’ procedure 535. For each of the candidate 
integer sets returned by the 'first baseline candidate integer 
set selection’ procedure 520, a candidate baseline is com- 
puted using least squares along with its associated least 
squares solution residual and processing enters an iterative 
block. A 'first baseline candidate quality checks’ procedure 
540 performs two quality checks on the particular candidate 
integer set. First is a baseline magnitude check which is 
satisfied when the baseline candidate computed by the 'least 
squares computation’ procedure 535 for a given candidate 
integer set is within a first predetermined range. In the 
illustrative embodiment, the first predetermined range is 
0.15 GPS LI wavelengths of the known magnitude of the 
baseline. Second is a least squares solution residual check 
that is satisfied when the magnitude of the least squares 
solution residual vector associated with the baseline candi- 
date falls within a second predetermined range. In the 
illustrative embodiment, the second predetermined range is 
0.13 GPS LI wavelengths. If either test is not satisfied, the 
baseline candidate is removed from further consideration. A 
'rank candidate’ procedure 545 then assigns a rank to 
baseline candidates passing both quality checks based on 
their least squares solution residual with higher rankings 
corresponding to smaller residuals. Iteration continues until 
a 'remaining candidates’ decision procedure 550 exits 
through its 'no’ branch indicating all of the baseline candi- 
dates have been considered. Processing continues to a 'pass- 
ing candidate’ decision procedure 555 that exits through its 
'no’ branch to execute a 'failure condition’ output block 575 
if no baseline candidate has passed the baseline quality 
checks. Otherwise the highest ranked and next highest 
ranked passing baseline candidates are returned to a 'differ- 
ence measurement’ procedure 560. The 'difference measure- 
ment’ procedure 560 computes a difference measure 
between the highest ranked baseline candidate and the next 
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highest ranked baseline candidate. In one embodiment of the 
invention the difference measure is the difference between 
the magnitudes of the least squares solution residual. A 
‘sufficient difference’ decision procedure 565 receives the 
difference measure. If the difference measure is above a 
predetermined threshold, the ‘sufficient difference’ decision 
procedure 565 exits through its ‘yes’ branch. In the illustra- 
tive embodiment, the predetermined threshold is 0.05 GPS 
LI wavelengths. When this occurs, a ‘first baseline and 
candidate integer set’ output block 570 returns the highest 
ranked baseline candidate and its associated candidate inte- 
ger set to the calling process 500 and processing completes 
through an ‘end’ terminal 580. If the difference measure is 
not above the predetermined threshold, the ‘sufficient dif- 
ference’ decision procedure 565 exits through its ‘no’ branch 
and the ‘failure condition’ output block 575 returns an 
indication to the calling process 500 that identification of the 
first baseline b 32 was unsuccessful based on the epoch of 
data currently under consideration. Processing then com- 
pletes through the ‘end’ terminal 580. 

The previous discussion of FIG. 5A, FIG. 5B, and FIG. 
5C described the backward integer search procedure to 
determine the first baseline in accordance with an illustrative 
embodiment of the invention. With reference to FIG. 4, 
processing has progressed until the ‘first baseline identified’ 
decision procedure 425 exits through its ‘yes’ branch. 

Because of the desire for larger antenna separations to 
produce more precise attitude estimates, in the illustrative 
embodiment, the fourth antenna A4 and the first antenna A1 
are located several wavelengths away from the second 
antenna A2 and the third antenna A3. This means that a 
single-epoch backward integer search procedure to resolve 
integer ambiguities associated with these antennas would be 
extraordinarily lengthy due to the vast number of integer set 
candidates, and also quite likely to result in incorrect integer 
set determination. Therefore, the successful conclusion of 
the ‘first baseline determination’ process 500, the integer 
search for the first baseline b 32 , is now capitalized upon to 
carry out a forward integer search procedure for the integers 
associated with the second baseline b 43 . A more complete 
description of the forward search for the second baseline 
follows. 

FIG. 6A illustrates a ‘second baseline determination’ 
process 600 invoked by the ‘forward search for second 
baseline’ procedure 430. Processing initiates at a ‘start’ 
terminal 610. Broadly speaking, the ‘second baseline deter- 
mination’ process 600 is similar to the ‘first baseline deter- 
mination’ process 500 discussed in reference to FIG. 5 A. 
The differences lie in the manner in which candidate integer 
sets are selected and the correct set identified by a ‘second 
baseline candidate integer set selection’ procedure 620 and 
a ‘second baseline correct candidate integer set identifica- 
tion’ procedure 630, respectively. A ‘successful identifica- 
tion’ decision procedure 685 exits through its ‘yes’ branch to 
a ‘second baseline and candidate integer set’ output block 
690 or through its ‘no’ branch to a ‘failure condition’ output 
block 695 and processing completes through an ‘end’ ter- 
minal 699 in an analogous manner to what was described 
with respect to the corresponding functional elements 
depicted in FIG. 5 A, the ‘successful identification’ decision 
procedure 585, the ‘first baseline and candidate integer set’ 
output block 590, and the ‘failure condition’ output block 
595. 

The ‘second baseline candidate integer set selection’ 
procedure 620 is shown in greater detail by FIG. 6B. 
Processing initiates at a ‘start’ terminal 621 and continues to 
a ‘first baseline’ input block 622 the presence of which 
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reflects the fact that the first baseline b 32 is used to aid the 
search for the second baseline b 43 . Next a ‘heading and roll 
estimation’ procedure 624 generates a coarse estimate of the 
heading and roll of the antenna configuration 200 based on 
5 the first baseline. More specifically, associated with the 
antenna configuration 200 is an antenna coordinate system 
consisting of a first coordinate unit vector Cl, a second 
coordinate unit C2 and a third coordinate unit vector C3. The 
second coordinate unit vector C2 is taken parallel to b 42 , the 
10 third coordinate unit vector C3 is taken parallel to the cross 
product b 42 xb 41 and the first coordinate unit vector is taken 
to complete a right-hand coordinate system. The ‘heading 
and roll estimation’ procedure 624 estimates heading and 
roll angles of the antenna coordinates with respect to local- 
15 vertical local-horizontal Earth-fixed coordinate system. One 
skilled in the art will recognize that selection of heading and 
roll from among the three orientation angles is not funda- 
mental to the present invention. 

Processing continues to a ‘second baseline limited geo- 
20 metric region determination’ procedure 626. Here, the head- 
ing and roll estimates are used in conjunction with the 
knowledge of the physical orientation arid length of the 
second baseline b43 to determine a limited geometric region 
for the location of as yet unknown endpoint of the second 
25 baseline b 43 , i.e. the location of the fourth antenna A4. The 
limited geometric region may be described as follows. 
Fixing the end of the second baseline located at the third 
antenna A3, rotate the baseline through a range of antenna 
configuration heading and roll values equal to the 
30 previously-computed heading and roll estimates plus/minus 
incremental offsets. The limited geometric region is the 
region swept out by the second baseline b 43 , i.e. a region on 
whose longitudinal axis lies an approximate location of the 
fourth antenna A4. Determination of the limited geometric 
35 region for use in the forward integer search procedure is 
further discussed in connection with FIG. 8 below. 

Next a ‘candidate selection from second baseline limited 
geometric region’ procedure 628 identifies candidate integer 
sets corresponding to points in the limited geometric region. 
40 In one embodiment of the invention the following relation- 
ships are used: 

DD^=[dLOS]b^/A, 

N^diD=round(DD i7s -DDj 

45 

where bars above a quantity indicates that it is a baseline 
candidate value, by is the second baseline b 43 appro- 
priately rotated by the previously-estimated antenna 
coordinate system heading and roll plus/minus an 
50 incremental offset and the term “round” indicates the 
quantity rounded to the nearest integer. By iterating 
through the range of offsets, candidate integer sets in 
the limited geometric region are identified and process- 
ing completes through an ‘end’ terminal 629. 

55 Depicted in FIG. 6C is the ‘second baseline correct 
candidate integer set selection’ procedure 630 which 
receives the candidate integer sets selected by the procedure 
620 previously described with reference to FIG. 6B. Pro- 
cessing initiates at a ‘start’ terminal 631 and continues to a 
60 ‘least squares computation’ procedure 635 that computes the 
least squares solution for the second baseline b 43 for can- 
didate integer sets from the limited geometric region. The 
least squares solution residual is also computed. Next, a 
‘second baseline candidate quality checks’ procedure 640 
65 performs slightly different quality checks on the second 
baseline candidates than those performed by the ‘first base- 
line candidate quality checks’ procedure 540. In addition to 
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the baseline magnitude check, and the least squares solution 
residual check, a score is computed for the second baseline 
candidate integer sets. The baseline magnitude check here is 
whether the magnitude of the second baseline b 43 is within 
0.11 GPS LI wavelengths of the known magnitude; the least 
squares solution residual check is wither the second baseline 
solution residual vector magnitude is less than 0.22 GPS LI 
wavelengths. Additionally, in the illustrative embodiment, 
the error in the calculated angle between the first baseline 
b 32 and the second baseline b 43 is checked to verify it is 
below 15.0 degrees, and the baseline score is checked to 
verify that it is below a value of 0.9. 

The score is a weighted sum of a baseline magnitude 
error, the least squares solution residual, and the absolute 
value of the angle error between the candidate the second 
baseline b 43 and the previously estimated first baseline b 32 . 
The baseline magnitude error is the difference (in GPS LI 
wavelengths) between the magnitude of the least squares 
solution for the second baseline candidate and the known 
magnitude of the second baseline b 43 . In the illustrative 
embodiment, the weights are 10, 1, and 0.1 for the least 
squares solution residual, baseline magnitude error, and 
absolute value of the angle error, respectively. 

One skilled in the art will recognize that imposition of the 
constraint that baseline candidates pass plural successive 
baseline quality tests improves the robustness of the baseline 
determination process. 

Execution of the remainder of the procedure 630 is similar 
to the corresponding procedure discussed in connection with 
FIG. 5C. Briefly, a 'rank candidate’ procedure 645 deter- 
mines the highest and next highest ranked candidate (based 
on the least squares solution residual) through iteration over 
the candidates controlled by a 'remaining candidates’ deci- 
sion procedure 650; a 'passing candidate’ decision proce- 
dure 655 reports a 'failure condition’ output block 675 if 
there are not passing candidates. Otherwise a 'difference 
measurement’ procedure 660 and a 'sufficient difference’ 
decision procedure 665 determine if the highest ranked 
candidate is sufficiently different from the next highest 
ranked candidate. In the illustrative embodiment, this dif- 
ference is whether the difference in the residual vectors of 
the best and second best candidates is greater than 0.05 GPS 
LI wavelengths. If so a 'second baseline and candidate 
integer set’ output block 670 reports the highest ranked 
candidate and its corresponding integer set; otherwise the 
'failure condition’ output block 675 is reported. Processing 
completes through an 'end’ terminal 680. 

Again making reference to FIG. 4, when processing has 
completed so the 'second baseline identified’ decision pro- 
cedure 435 exits through its 'yes’ branch, the 'forward 
search for third baseline’ procedure 440 initiates. Execution 
of the procedure 440 is similar to the 'forward search for 
second baseline’ procedure 430 except for the determination 
of the candidate integer sets. 

FIG. 7A illustrates the 'third baseline determination’ 
process 700 in greater detail. Processing initiates at a 'start’ 
terminal 710 and continues to a 'third baseline candidate 
integer set selection’ procedure 720 and a 'third baseline 
correct candidate integer set identification’ procedure 730 
that will be discussed in connection with FIG. 7B and FIG. 
7C, respectively. Execution of the remainder of the 'third 
baseline determination’ process 700 is similar to that of the 
procedures discussed previously in connection with FIG. 5 A 
and FIG. 6A. Briefly, a 'successful identification’ decision 
procedure 785 tests whether the third baseline b 21 was 
successfully identified. If so a 'third baseline and candidate 
integer set’ output block 790 reports the third baseline b 21 
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and the corresponding integer set; otherwise a 'failure 
condition’ output block 795 executes for this epoch of data. 
Processing completes through an 'end’ terminal 799. 

Referring now to FIG. 7B, processing initiates at a 'start’ 
5 terminal 721 and a 'second baseline’ input block 722 
receives the second baseline b 43 determined as described 
above in connection with FIG. 6 A, FIG. 6B, and FIG. 6C. A 
'heading and roll estimation’ procedure 724 uses the second 
baseline b 43 to re-estimate the heading and roll angles of the 
10 antenna coordinate with respect local-vertical and local- 
horizontal Earth-fixed coordinates. Processing continues to 
a 'third baseline limited geometric region determination’ 
procedure 726 and a 'candidate selection from third baseline 
limited geometric region’ procedure 728. These procedures 
15 execute in a manner analogous to the corresponding proce- 
dures discussed in connection with FIG. 6B with the fol- 
lowing exception: the third baseline limited geometric 
region is nearly planar. It corresponds to the region swept 
out by the third baseline b 21 adjacent to the second antenna 
20 A2 as antenna configuration pitch angles are varied through 
an appropriate angle range, for example +/-20 degrees. The 
width of the nearly planar region searched may be zero or 
may encompass a width of several degrees, constructed by 
incremental offsets in heading. Candidate integer sets are 
25 chosen from within this region in a manner similar that 
described in connection with FIG. 6B except the pitch and 
heading angles are offset rather than the heading and roll 
angles. One skilled in the art will recognize the concept 
described above is not limited to pitch and heading, but also 
30 encompasses other orientation angles. Processing completes 
through an 'end’ terminal 729. 

The third baseline candidate integer sets having been 
selected, the 'third baseline determination’ process 700 
continues to a 'third baseline correct candidate integer set 
35 identification’ procedure 730 as depicted in FIG. 7C. This 
procedure executes analogously to the 'second baseline 
correct candidate integer set identification’ procedure 630 
and reference should be made to the discussion of the 
procedure 630 for a description of the 'third baseline correct 
40 candidate integer set identification’ procedure 730. Briefly, 
processing initiates at a 'start’ terminal 731 and continues to 
a 'least squares computation’ procedure 735 for computing 
the least squares solution and residuals of baseline candi- 
dates. An iterative block controlled by a 'remaining candi- 
45 dates’ decision procedure 750 executes baseline quality tests 
at a 'third baseline candidate quality checks’ procedure 740 
and a 'rank candidate’ procedure 745 assigns rank to can- 
didates on the basis of a ranking criterion, for example the 
magnitude of the least squares residual. In the illustrative 
50 embodiment, the third baseline quality checks are whether 
the difference in magnitude of the third baseline from the 
known magnitude is less than 0.11 GPS LI wavelengths; 
whether the magnitude of the residual vector is less than 
0.10 GPS LI wavelengths; whether the score for the third 
55 baseline is less than 1.3; and whether the error in the angle 
calculated between the second and third baselines is less 
than 7.0 degrees. When iteration is complete, a 'failure 
condition’ output block 775 executes for this epoch of data 
if a 'passing candidate’ decision procedure 755 determines 
60 no candidates passed the baseline quality tests. Otherwise a 
'difference measurement’ procedure 760 and a 'sufficient 
difference’ decision procedure 765 function to test whether 
a highest ranked candidate and a next highest ranked can- 
didate are sufficiently different. In the illustrative 
65 embodiment, the residual difference is tested to verify it is 
above 0.10 GPS LI wavelengths. If so, the 'sufficient 
difference’ decision procedure 765 exits through its 'yes’ 
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branch and a 'third baseline and candidate integer set’ output 
block 770 reports the third baseline b 21 and the correspond- 
ing integer set; otherwise the 'failure condition’ output block 
775 executes for this epoch of data. Processing completes 
through an 'end* terminal 780. 

FIG. 8 illustrates a plurality of baseline candidate vectors 
in a limited geometric region 800 in accordance with the 
forward integer search procedures previously described in 
reference to FIG. 6B and FIG. 7B. The plurality of baseline 
candidate vectors in a limited geometric region 800 further 
illustrates the nature of the forward integer search procedure. 
As depicted, a forward integer search is carried out to 
determine the baseline between a first antenna A,- and a 
second antenna Ay. A correct baseline b~ is depicted. The 
correct baseline b tj of one of a plurality of candidate base- 
lines b-j. While the limited geometric region is depicted as 
having circular cross sections in planes perpendicular to a 
line joining the first antenna A,- and the second antenna Ay, 
one skilled in the art will recognize this is not fundamental 
to the current invention. Embodiments of the invention 
alternatively use regions having rectangular cross-sections. 

The forward integer search procedure uses information 
respecting one or more orientation angles of the baseline 
being determined. This information may come from previ- 
ously determined baselines or from other sources. In accor- 
dance with the illustrative embodiment, the orientation angle 
information are coarse estimates of the heading and roll of 
the correct baseline b tj from one or more previously deter- 
mined baselines. The limited geometric region is selected 
responsive to the orientation angle information and the 
known magnitude of the correct baseline b iy . 

Shown in the limited geometric region are the plurality of 
candidate baselines b ( y. The plurality of candidate baselines 
b-j correspond to candidate integer sets of the integer ambi- 
guity associated with the baseline between the first antenna 
A,- and the second antenna Ay. One skilled in the art will 
recognize that not all candidate integer sets would have 
corresponding baseline candidates in the limited geometric 
region. Baseline candidate vectors lying in the limited 
geometric region are evaluated by the forward integer search 
procedure. The forward integer search procedure thus 
selects the correct baseline b^ from among the plurality of 
candidate baselines b tj lying in the limited geometric region. 

EXAMPLES 

Any single baseline has a pointing direction, including 
any baseline derived from any combination of two or more 
antennas. This pointing direction can be used to determine 
partial orientation of a rigid body, by calculation of two or 
three degrees of freedom of the rigid body. One example is 
the determination of heading and roll of a rigid body to 
which the two or more antennas are mounted. Another 
example is determination of heading and pitch. The partial 
determination of orientation is not limited to these examples. 

Any non-parallel set of baselines, comprising two or more 
baselines, derived from three or more antennas can be used 
to completely determine the orientation of a rigid body to 
which the antennas are attached. Orientation can be 
expressed in a number of ways, one of which is by use of 
orientation angles such as heading, pitch, and roll. 

Partial orientation or complete orientation of a rigid body 
can be used for information purposes and for operations 
including control and navigation of vehicles, aircraft, water- 
craft and spacecraft. 

In similar fashion, partial or complete orientation infor- 
mation can be used for determination of pointing direction 
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for sighting devices and remote sensing devices. Orientation 
information can be used for pointing devices intended to 
launch projectiles, both ballistic and propelled. 

Pointing direction and orientation information can be used 
5 for identifying the physical shape and location of structures, 
rigid or non-rigid, to which multiple antennas are attached. 
This is useful, for example, in determining the shape and 
location of portions of flexible or rigid surfaces or booms, or 
that of implements used to carry out physical placement or 
10 excavation. This information can be used among other uses, 
for controlling the shape and position of such flexible or 
rigid surfaces, booms, or implements. 

Orientation and pointing information can also be used to 
aim devices designed to receive electromagnetic radiation, 
15 such as antennas and solar panels. 

Orientation and pointing information can be used to 
determine the surface characteristics of a surface over which 
an instrument is moving. 

20 While the present invention may be used alone as the 
basis for an attitude determination system, it may be com- 
bined with other known like methods, such as inertial 
sensors. 

Although the present invention has been described in 
25 terms of an illustrative embodiment, one skilled in the art 
will understand that various modifications and alterations 
may be made without departing from the scope of the 
invention. Accordingly, the scope of the invention is not to 
be limited to the particular embodiments discussed herein, 
30 but should be defined only by the appended claims and 
equivalents thereof. 

What is claimed is: 

1. A method for determining inter- antenna baselines from 
radio wave phase measurements, said method comprising: 
35 identifying a first antenna and a second antenna located a 
distance D1 apart; 

identifying a third antenna located a distance D2 from the 
second antenna; 

determining a first baseline between the first antenna and 
40 the second antenna using a backward integer search 
procedure; and 

determining a second baseline between the second 
antenna and the third antenna using a forward integer 
search procedure responsive to the first baseline. 

45 2. The method of claim 1, wherein the distance D1 is less 

than or equal to 2.5 wavelengths of a carrier wave that 
emanates from a radio source. 

3. The method of claim 1, wherein the distance D2 is 
greater than or equal to 1.0 wavelengths of a carrier wave 

50 that emanates from a radio source. 

4. The method of claim 1, further comprising the steps: 

establishing a fourth antenna; and 

determining a third baseline between the fourth antenna 
55 and one of the first antenna, the second antenna, and the 
third antenna, the third baseline being determined using 
a forward integer search procedure responsive to at 
least one previously-determined baseline. 

5. The method of claim 1, wherein the step of determining 
60 a first baseline between the first antenna and the second 

antenna using a backward integer search procedure com- 
prises steps of: 

performing a backward integer search procedure for each 
of a plurality of different radio source carrier frequen- 
65 cies; 

identifying at least one baseline candidate at each of the 
different radio source carrier frequencies; 
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evaluating each of the baseline candidates with a baseline 
quality test; and 

selecting as a first baseline the baseline candidate that 
satisfies the baseline quality test for each of the plu- 
rality of different radio source carrier frequencies. 

6. The method of claim 1, wherein the step of determining 
a second baseline between the second antenna and the third 
antenna using a forward integer search procedure responsive 
to the first baseline comprises steps of: 

performing a forward integer search procedure for each of 
a plurality of different radio source carrier frequencies; 

identifying at least one baseline candidate at each of the 
different radio source carrier frequencies; 

evaluating each of the baseline candidates with a baseline 
quality test; and 

selecting as a first baseline the baseline candidate that 
satisfies the baseline quality test for each of the plu- 
rality of different radio source carrier frequencies. 

7. The method of claim 1, wherein the step of determining 
a first baseline between the first antenna and the second 
antenna using a backward integer search procedure com- 
prises steps of: 

ranking a plurality of baseline candidates on the basis of 
a ranking criterion; 

identifying a highest ranked baseline candidate and a next 
highest ranked baseline; 

computing a difference measure related to the difference 
between the highest ranked baseline candidate and the 
next highest ranked baseline candidate; and 

selecting the highest ranked baseline candidate as the first 
baseline if the difference measure is above a predeter- 
mined threshold. 

8. The method of claim 1, wherein the step of determining 
a second baseline between the second antenna and the third 
antenna using a forward integer search procedure responsive 
to the first baseline comprises steps of: 

determining a plurality of attitude component estimates 
from the first baseline; 

selecting a limited geometric region responsive to the 
plurality of attitude component estimates, vectors in the 
limited geometric region corresponding to candidate 
integer sets; 

determining a second baseline responsive to the candidate 
integer sets corresponding to vectors in the limited 
geometric region. 

9. A method for selecting a baseline from a plurality of 
baseline candidates comprising: 

ranking the plurality of baseline candidates on the basis of 
a ranking criterion; 

identifying a highest ranked baseline candidate and a next 
highest ranked baseline; 

computing a difference measure related to the difference 
between the highest ranked baseline candidate and the 
next highest ranked baseline candidate; and 

selecting the highest ranked baseline candidate if the 
difference measure is above a predetermined threshold. 

10. The method of claim 9, wherein the difference mea- 
sure is the difference of the magnitude of the least squares 
residuals of the baselines associated with the highest ranked 
candidate integer set and the next highest ranked candidate 
integer set. 

11. A method of determining an inter- antenna baseline 
from estimates of two or more orientation angles of the 
baseline, the method comprising: 
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determining a limited geometric region from estimates of 
the two or more orientation angles of the baseline; 
identifying a plurality of baseline candidate vectors in the 
limited geometric region from candidate integer sets 
? that resolve an integer ambiguity associated with the 
inter-antenna baseline; 

selecting the inter-antenna baseline from the plurality of 
baseline candidate vectors. 

10 12. The method of claim 11, wherein the limited geomet- 

ric region is determined by deviating the estimates of the two 
or more orientation angles by a predetermined amount. 

13. The method of claim 11, wherein the step of selecting 
the inter-antenna baseline from the plurality of baseline 

15 candidate vectors comprises steps of: 

ranking the plurality of baseline candidates on the basis of 
a ranking criterion; 

determining a highest ranked candidate and a next highest 
2Q ranked baseline candidate; 

computing a difference measure related to the difference 
between the highest ranked baseline candidate and the 
next highest ranked baseline candidate; and 
rejecting the plurality of baseline candidates if the dififer- 
25 ence measure is within a predetermined range; and 
selecting the inter-antenna baseline from the plurality of 
baseline candidates if difference measure is outside the 
predetermined range. 

14. The method of claim 13, wherein the difference 
30 measure is the difference of the magnitude of the least 

squares residuals of the baseline candidates. 

15. The method of claim 11, wherein the step of selecting 
the inter-antenna baseline from the plurality of baseline 
candidate vectors comprises steps of: 

or 

evaluating the plurality of baseline candidate vectors with 
a baseline quality test; and 

selecting the inter-antenna baseline from baseline candi- 
date vectors passing the baseline quality test. 

40 16. The method of claim 15, wherein the baseline quality 

test is based on the baseline’s magnitude. 

17. The method of claim 15, wherein the baseline quality 
test is based on the baseline’s least squares residual. 

18. The method of claim 15, wherein the baseline quality 
45 test is based on a weighted sum of the errors in the baseline’s 

computed magnitude, the magnitude of the baseline’s least 
squares residual, and an error in the computed angle formed 
with the first baseline, wherein some weights may be zero. 

19. An apparatus configured to determine inter- antenna 
5Q baselines from radio wave phase measurements; the appa- 
ratus having a processor, and a memory, the apparatus 
comprising: 

means for determining a first baseline between a first 
antenna and a second antenna using a backward integer 
55 search procedure, wherein the first antenna and the 
second antenna are located a distance D1 apart; and 
means for determining a second baseline between the 
second antenna and a third antenna using a forward 
integer search procedure responsive to the first 
60 baseline, wherein the second antenna and the third 
antenna are located a distance D2 apart. 

20. The apparatus of claim 19, wherein the distance D1 is 
less than about 2.5 wavelengths of a carrier wave that 
emanates from a Global Positioning System satellite. 

65 21. The apparatus of claim 19, wherein the distance D2 is 

greater than about 1.0 wavelengths of a carrier wave that 
emanates from a pseudolite. 
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22. The apparatus of claim 19 further comprising: 
means for determining a third baseline between a fourth 

antenna and one of the first antenna, the second 
antenna, and the third antenna, the third baseline being 
determined using a forward integer search procedure 5 
responsive to at least one previously-determined base- 
line. 

23. The apparatus of claim 19 wherein the means for 
determining a first baseline between the first antenna and the 
second antenna using a backward integer search procedure 10 
comprises: 

means for ranking a plurality of baseline candidates on the 
basis of a ranking criterion; 

means for identifying a highest ranked baseline candidate 
and a next highest ranked baseline; 
means for computing a difference measure related to the 
difference between the highest ranked baseline candi- 
date and the next highest ranked baseline candidate; 
and; 2 o 

means for selecting the highest ranked baseline candidate 
as the first baseline if the difference measure is above 
a predetermined threshold. 

24. The apparatus of claim 19, wherein the means for 
determining a second baseline between the second antenna 25 
and the third antenna using a forward integer search proce- 
dure responsive to the first baseline comprises: 

means for determining a plurality of attitude component 
estimates from the first baseline; 
means for selecting a limited geometric region responsive 30 
to the plurality of attitude component estimates, vectors 
in the limited geometric region corresponding to can- 
didate integer sets; 

means for determining a second baseline responsive to the 35 
candidate integer sets corresponding to vectors in the 
limited geometric region. 

25. The apparatus of claim 23 wherein the difference 
measure is the difference of the least squares residuals of the 
baselines associated with the highest ranked candidate inte- 4Q 
ger sets and the next highest ranked candidate integer sets. 

26. An apparatus for determining inter-antenna baselines 

among at least a first antenna, a second antenna, and a third 
antenna from a plurality of carrier waves, the apparatus 
comprising: 45 

a microprocessor; 

at least one receiver configured to report the carrier wave 
phases sensed by the antennas; 
a storage medium coupled to the microprocessor and 
storing code which configure the microprocessor to 50 
perform steps of: 

detecting a sensed phase difference between the phases 
reported by the first receiver and the second receiver for 
each of the plural carrier waves; 
identifying candidate integer sets constituting an integer 
ambiguity associated with a first baseline, the first 
baseline representing the relative position of the first 
antenna and the second antenna; 
determining the first baseline from the sensed phases 60 
differences and the candidate integer sets wherein the 
integer ambiguity is resolved with a backward integer 
search procedure; and 

determining a second baseline representing the relative 
position of the second antenna and the third antenna; 65 
using a forward integer search procedure responsive to 
the first baseline. 
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27. An apparatus for selecting a baseline from a plurality 
of baseline candidates, the apparatus comprising: 

a microprocessor; 

a storage medium coupled to the microprocessor and 
storing program instructions that configure the micro- 
processor to perform steps of: 

ranking the plurality of baseline candidates on the basis 
of a ranking criterion; 

identifying a highest ranked baseline candidate and a 
next highest ranked baseline; 
computing a difference measure related to the differ- 
ence between the highest ranked baseline candidate 
and the next highest ranked baseline candidate; and 
selecting the highest ranked baseline candidate if the 
difference measure is above a predetermined thresh- 
old. 

28. An apparatus for determining an inter- antenna base- 
line from estimates of two or more orientation angles of the 
baseline, the apparatus comprising: 

a microprocessor; 

a storage medium coupled to the microprocessor and 
storing program instructions that configure the micro- 
processor to perform steps of: 

determining a limited geometric region from estimates 
of the two or more orientation angles of the baseline; 
identifying a plurality of baseline candidate vectors in 
the limited geometric region from candidate integer 
sets that resolve an integer ambiguity associated with 
the inter- antenna baseline; 

selecting the inter-antenna baseline from the plurality 
of baseline candidate vectors. 

29. The apparatus of claim 28, wherein the step of 
selecting the inter-antenna baseline from the plurality of 
baseline candidate vectors comprises steps of: 

ranking the plurality of baseline candidates on the basis of 
a ranking criterion; 

identifying a highest ranked baseline candidate and a next 
highest ranked baseline; 

computing a difference measure related to the difference 
between the highest ranked baseline candidate and the 
next highest ranked baseline candidate; and 
selecting the highest ranked baseline candidate if the 
difference measure is above a predetermined threshold. 

30. The apparatus of claim 29, wherein the step of 
selecting the inter-antenna baseline from the plurality of 
baseline candidate vectors comprises steps of: 

evaluating the plurality of baseline candidate vectors with 
a baseline quality test; and 

selecting the inter-antenna baseline from baseline candi- 
date vectors passing the baseline quality test. 

31. The apparatus of claim 30, wherein the baseline 
quality test is based on the baseline’s magnitude. 

32. The apparatus of claim 30, wherein the baseline 
quality test is based on the baseline’s least squares residual. 

33. The apparatus of claim 30, wherein the baseline 
quality test is based on a weighted sum of the errors in the 
baseline’s computed magnitude, the magnitude of the base- 
line’s least squares residual, and an error in the computed 
angle formed with the first baseline, wherein some weights 
may be zero. 

34. An antenna configuration comprising 
a first antenna 

a second antenna located a distance D1 from the first 
antenna, wherein the distance D1 is less than or equal 
to 1.75 wavelengths of a predetermined carrier wave; 
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a third antenna located a distance D2 from the second 
antenna, wherein the distance D2 is greater than or 
equal to 2.75 wavelengths of the predetermined carrier 
wave; and 

a fourth antenna located a distance D3 from said first 
antenna and a distance D4 from said third antenna, 
wherein said distance D3 is greater than or equal to 3.0 
wavelengths of the predetermined carrier wave and 
wherein said distance D4 is greater than or equal to 1.0 
wavelengths of the predetermined carrier wave. 

35. The antenna configuration of claim 34 wherein the 
predetermined carrier wave is a carrier wave from a Global 
Positioning System satellite. 

36. A GPS attitude system comprising 

a plurality of antenna-receiver units configured to report 
phases from GPS carrier waves; and 

a processing system coupled to the plurality of antenna- 
receiver units configured to process the phases to 
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compute attitude, the processing system configured to 

initialize with steps comprising: 

identifying a first antenna, a second antenna, and a third 

antenna associated with the plurality of antenna- 
5 • 

receiver units; 

determining a first baseline between the first antenna 
and the second antenna using a backward integer 
search procedure; and 

io determining a second baseline between the second 
antenna and the third antenna or between the first 
antenna and the third antenna using a forward integer 
search procedure responsive to the first baseline. 

37. The antenna configuration of claim 34 wherein said 
fourth antenna is noncolinear with said first antenna and said 
third antenna. 



